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Abstract. We analyze the ground state of a strongly interacting fermion chain with a 
supersymmetry. We conjecture a number of exact results, such as a hidden duality between weak 
and strong couplings. By exploiting a scale free property of the perturbative expansions, we find 
exact expressions for the order parameters, yielding the critical exponents. We show that the ground 
state of this fermion chain and another model in the same universality class, the XYZ chain along a 
line of couplings, are both written in terms of the same polynomials. We demonstrate this explicitly 
for up to N = 24 sites, and provide consistency checks for large N. These polynomials satisfy a 
recursion relation related to the Painleve VI differential equation, and using a scale-free property of 
these polynomials, we derive a simple and exact formula for their N — ¥ oo limit. 



PACS numbers: 11.30.Pb, 04.60.Nc 

1. Introduction 

The calculation of exact results for physical quantities in interacting one-dimensional quantum chains 
has a long and distinguished history, going back to the early days of quantum mechanics. Bethe's 
results for the Heisenberg chain more than 75 years ago were the first; there are now a whole host of 
models that are often referred to as "solved" . 

The word "solved" sometimes means that the exact ground state is known in closed form, and it 
sometimes means that the model is integrable. Integrability means that there exist enough conserved 
charges to completely constrain the dynamics. Although it is often possible to exploit the integrability 
to do exact computations, usually it is not possible to find the ground state explicitly in an integrable 
model. Conversely, many models where the ground state is known explicitly are not integrable. 

Even when we restrict our study to integrable models, there is still a wide variation in what can 
be calculated and how difficult it is to do so. On one end of the spectrum, there are systems that can 
be mapped to free fermions, and so essentially everything can be computed explicitly. The canonical 
example of this is the quantum Ising chain, often referred to as the Ising model in a transverse field. 
On the other end of the spectrum, there are models such as the XYZ chain, where the standard Bethe 
ansatz does not work, and much more elaborate methods are required to derive exact results pQ. 

It is our aim in this paper to discuss one example of a class of integrable models that are seemingly 
difficult to solve, in that the standard Bethe ansatz does not apply. In fact, in the model discussed 
in this paper, we have not even proved that it is integrable (although we believe it is), much less 
exploited the integrability directly. However, we will show that exact results can be found by using 
very different methods than those usually utilized in integrability. 

An essential ingredient in our analysis is supersymmetry. One of the hallmarks of supersymmetric 
systems is the special properties satisfied by the ground state [2]. This usually makes it possible 
to not only find the ground-state energy exactly and easily for any number of sites N, but to find 
simple operators that annihilate the ground state. A supersymmetric system need not be integrable, 
but when it is, the ground state itself satisfies a number of special properties, even for an integrable 
system. For example, we will show here that this makes it possible to find exact values for the critical 
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exponents simply by analyzing the ground state at finite N. It also makes it possible to numerically 
determine scaling functions extremely accurately from very small system sizes. 

The specific systems we are studying were introduced in [3] . These are chains of spinless fermions 
with strong interactions where the supersymmetry is built in from the start, so that the ground state of 
this many-body system has many special properties. The supersymmetric fermion chains analyzed in 
depth in [3j |4] are critical. Here we study a one-parameter deformation of the critical supersymmetric 
chain. As we discussed briefly in an earlier paper [5] , we find simple but exact expressions for a variety 
of physical quantities. Here we extend this work considerably, and find both scaling functions and 
exact results for the ground state itself. 

A famous example of a system in the class we are studying is the XXZ chain at anisotropy 
A = —1/2 (see e.g. [6j [3 [8]). A host of elegant results relating the ground state of this system to 
combinatorial quantities have been found. One of the reason these properties were found is that like 
with a supersymmetric system, the ground state energy is exactly known. In fact, this XXZ chain 
has a hidden supersymmetry [H [3] . Our off-critical results on the fermion chain here results are quite 
analogous to those found for an off-critical deformation of the XXZ chain to the XYZ chain studied 
in [TU1 [TTJ [12l [T3l E] . In fact, not only is this XYZ chain in the same universality class as our fermion 
chain, but the ground states can be constructed in terms of the same polynomials! These polynomials 
have a host of very special properties that we analyze in depth. 

We believe that supersymmetric integrable systems of this type are "more than integrable" - 
the additional symmetries coming from supersymmetry result in structure simply not present in an 
integrable system. In fact, in some aspects, supersymmetric integrable systems are simpler than free- 
fermion systems. For example, the ground state energy in a free-fermi system such as the Ising chain 
depends non-trivially on the size of the system: as one increases N, one must keep filling the fermi 
sea. In our supersymmetric systems, the ground state energy is exactly zero for any N. 

The layout of this paper is as follows. In section [5] we introduce the lattice fermion model with 
supersymmetry. Using both the supersymmetry and an exact computation of the ground state up to 
N = 24, we discuss general and specific properties of the ground state in section [3j In particular, we 
show how the ground states can be characterized by a polynomial with very special properties. In 
section |4j we discuss the duality and the scaling behavior of the order parameters. A connection to 
the XYZ model, discussed in section [5j indicates a relation of our model with classically integrable 
equations. We exploit this to find a recursion relation for the polynomials describing the ground state. 
We then are able to find new results for these polynomials by exploiting their scale-free behavior. We 
present our conclusions along with some open questions in section [6] 

2. The model 

In this section we introduce our model explicitly, and analyze some of its basic properties. The 
degrees of freedom are "fat" spinless fermions on a periodic chain. They are annihilated and created 
by operators a and c] obeying the standard anticommutation relations {a,Cj} — {cj,cj} = and 
{ci,Cj} = Sij. What makes them fat is that, following [3], we do not allow two adjacent sites to be 

occupied. Defining the projector P 3 ■ = 1 — CjCj, this means we restrict the Hilbert space to include 
only configurations where PjPj + \ = for all j. Denoting an empty site by o, and an occupied site by 
• , this hard-core constraint therefore excludes pairs 

The reason we study fat fermions is that there exists a Hamiltonian with J\f = 2 supersymmetry 
on this Hilbert space [3]. This is because the "supercharges" 

N N 

are both nilpotent: (Q~) 2 = (Q + ) 2 = for any complex numbers Xj. Here we take periodic boundary 
conditions on a chain of length N, so all indices are interpreted mod N. Defining the Hamiltonian to 
be the anticommutator H = {Q + , the nilpotency requires that both Q~ and Q + commute with 
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H. Moreover, the fermion anticommutators ensure that H is local: 

N 

H = H [^-i (ViA*c] +lCj + h.c.) P j+2 + \X 3 \ 2 P ^P J+l ] (1) 

i=i 

The first term is thus a hopping term that preserves the hard-core exclusion. The second term can be 
interpreted as a potential, yielding a contribution |Aj| 2 whenever the sites j—1 andj + 1 are both empty. 
Since Q ± change the fermion number by ±1, H preserves the total number of fermions. In terms of 
operators, the fermion number operator F = Ylf=i c ] c j obeys [Q^-F 1 ] = ±Q ± and [Q ± ,H] = 0. 
The ensemble of algebraic relations between the four operators H, Q ± , F is the well-known N = 2 
supersymmetry algebra [2]. 

The supersymmetry has considerable consequences for the spectrum of H . Some important ones 
are: (i) all energies are positive or zero E > 0, (ii) all eigenstates with E > can be organized into 
doublets of the supersymmetry algebra, (Hi) an eigenstate has E = if and only if it is annihilated 
by both Q + and Q~ . Because of the latter property, finding the exact number of E = ground states 
is equivalent to computing the cohomology of Q, which is the dimension of the space of states that 
are annihilated by Q but cannot be written as Q acting on something else [2]. This computation is 
simple to do in this one-dimensional case [3] . The Hamiltonian (JlJ with periodic boundary conditions 
has exactly two zero-energy ground states when the number of sites is a multiple of three and all the 
Xj obey < |Aj| < oo. Moreover, the number of fermions in both these ground states is n = N/3. 
Henceforth we restrict our analysis to N = 3n sites and study these ground states in depth. 

The case Xj — 1 for all sites j was extensively studied in [3l 0] . The model is solvable by means of 
the coordinate Bethe ansatz. The Bethe equations are similar to those for the XXZ model at A = —1/2 
with twisted boundary conditions [5J [5] , and indeed in given momentum sector a mapping between 
the two models can be constructed [H The standard Bethe ansatz analysis of the excitations 
above these ground states shows that the spectrum is gapless, and so the model is critical in this 
case. The field theory discribing the scaling limit is a simple free massless boson, compactified at the 
supersymmetric radius and therefore the simplest field theory with J\f = (2, 2) supersymmetry (in the 
continuum the supersymmetry is enhanced to two chiral and two anti-chiral supersymmetries). 

When all Xj = 1, the model possesses a translation symmetry. Defining the translation operator 
T by T _1 c-T = ct +1 , we have [H,T] — 0. The cohomology computation shows that translation 
invariance is spontaneously broken when the number of sites N — in is a multiple of three: the 
eigenvalues of T on the two ground states are (— l) n e ±l7T ^ 3 . Thus translation by three sites T 3 on the 
lattice turns into translation symmetry in the field theory describing the scaling limit. (In the field 
theory, the eigenvalues of T in general are related to those of the Z3 subgroup of the chiral fermion 
number Q3].) 

Staggering the couplings typically destroys criticality in lattice chains. Here this is the case as 
well. However, since the model remains supersymmetric for any choice of the Xj, it is possible to 
preserve supersymmetry in the non-critical theory. We choose staggered coupling constants Xj with a 
period of three sites Xj + ^ p = Xj, in order to keep the same translation symmetry in the scaling limit. 
To preserve parity symmetry, i.e. invariance under reversal of the order of sites, we choose 

= ^3i+2 = V, A3i + i = 1 . (2) 

Defining PcjP^ 1 = cx-j+i, we have [H, P] = and [T 3 ,H] = 0. For simplicity, we focus on the 
case of real coupling constants y. Then it follows that in the staggered case the supercharges are 
linear functions of the coupling constant ~ + yQf, and the Hamiltonian becomes a quadratic 
polynomial H = Hq + yPi \ + y 2 H2- 

There is only one Lorentz-invariant supersymmetry-preserving deformation of the free boson 
conformal field theory, the sine-Gordon field theory at the supersymmetric point (/3 2 — 167r/3 in 
conventional normalization; for an overview see |15jh We thus expect the scaling limit of this staggered 
lattice model to be described by this massive field theory. This field theory is integrable, and we expect 
that the quantum chain is as well. This assertion is easy to check by finding the spectrum of H as 
a function of y on the computer for finite sizes. One finds that the characteristic polynomial of H 
factorizes into polynomials in y 2 with integer coefficients, characteristic of an enhanced symmetry. 
Moreover, there are many level crossings as y is varied, again characteristic of integrability. Although 
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we will not exploit this integrability directly in our analysis, the special properties we find for this 
model certainly stem from its presence. 



3. Properties of the ground states 

Because we know the ground-state energy is exactly zero for any N, it is natural to expect that the 
ground state itself will have special properties. In a supersymmetric theory one expects even more, 
since Q + and Q~ individually annihilate the ground state. In this section we describe first how this 
automatically requires some useful properties of the ground states. We then show how these ground 
states possess a number of much more surprising and elegant properties. 

The supersymmetry means it is easy to find the quantum numbers of the ground states for 
any value of the staggering. The cohomology computation shows that the two ground states have 
eigenvalue of (— l) N+1 under T 3 , while one is parity even and the other parity odd. We thus label the 
two ground states as \^ c (y)) and \$>°(y)). When the superscript is omitted, the equation applies to 
both ground states. 



3.1. Asymptotic limits 

One nice feature of introducing the staggering is that the model can easily be solved when y — » and 
y — > oo. This provides very useful intuition into the general situation, and an example of the power 
of supersymmetry. 

In the limit y — > 0, the ground state must solve the equations Qq |W) = Qq I*) = 0- For Qq to 
annihilate \ty) it must not contain any particles on the sites 2,5,8, .. . . For Qq to annihilates the 
ground state as well, any such site must have a particle to its left or right; the hard-core constraint then 
forbids creating a particle on these sites. There are only two configurations with these two properties, 
each with n = N/3 fermions: 

\a) =|oo.oo....)= m c u io>, 



P\A) = \* 00.00...)= (JI C 1- 2 J |0>- 

where |0) is the empty state, and to eliminate an annoying sign we define P|0) = (— l) Ar+1 |0). Parity 
eigenstates are given by the linear combinations 

\^(y = 0)) = \A)+P\A), \*°(y = 0)) = \A)-P\A) . (3) 

In the limit y — > +00, the ground state must be a solution of Qi\^) = Q^\^) = 0. We know that 
both Qi and Qf annihilate and create fermions only on sites 1, 3, 4, 6, . . . If all the sites 2,5,8,... are 
occupied, the hard-core constraint means it is impossible to create or annihilate particles on any other 
sites. Thus the polarized state 

\* e (y -+ 00)) a |P) = I o . o o . o • • •) = (f[ 4 j |0) 

is a ground state with even parity in this limit. In order to construct the parity-odd ground state, we 
notice that any state with n particles and all sites 2,5,8,... empty contains exactly one particle on 
each pair of sites 3j, 3j + 1 . It is therefore annihilated by Qf . In order for the state to be annihilated 
by Qj~ the particles on these two sites must occur as "singlets" | ■■■ o •■■■) — | o ■■ ■), and thus the 
second ground state appears as tensor product of such singlets 



n 

-+ co) cx |S) = [] (4j+i 4 3 ) 10) • 
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3.2. Characterizing by the ground state by a polynomial 

Since the ground states are annihilated by H for all y, we can imagine constructing the ground state 
by iteration. Namely, H = H Q + yHi + y 2 H 2 , where H = {Qq,Qq}. We just saw that |^(0)) is 
annihilated by both Qq and Q$ , so that Ho\^(0)) = 0. Then we define 

\9(y)) = |* )+y|* 1 )+y 2 |* 2 ) + ... (4) 

where the Vt^ are independent of y. Since H\^>(y)) = 0, we therefore have 

flo|*i) = -fli|*o) 

H \* a ) = - ffi|*«_i> - H 2 \* a _ 2 ) . (5) 

for a > 2. Thus if we could invert Hq, then by iterating we can determine any from |Wo). 

Of course, Hq is not invertible on the full Hilbert space, since Hq\^o) = 0. However, it is easy to 
check that the states given in © are the unique states annihilated by Hq for a given parity. Thus let 
us write for the N = 3n-site system, 

\9(y)) = m(y)\9 ) + \9), 

where ( x I'o|^ / ) = 0. Then given m(y), this iteration procedure uniquely determines \^S/). Moreover, 
since Ho is diagonal, it is easily inverted in the Hilbert space with jvPo) excluded. 

The ground state can thus be characterized by this function m(y). We will conjecture a recursion 
relation for it in n by noting the connection to the XYZ chain and the Painlcvc VI differential equation. 
Because this analysis is technically somewhat involved, we defer this calculation to section[5] However, 
several important properties of m(y) and the ground states follow solely from the supersymmetry. 

To illuminate these properties, it is useful to define the operator 



3jH 



L c 3j + 1 (6) 



3=1 

that counts the number of fcrmions on the "staggered" sites 2,5,8,... Since Qf does not change 
this number, [S, Qf ] = 0, but [S, Q^} = ±Q^- Therefore S commutes with Hq and H 2 , but not 
H\. The latter hops a particle on or off one of the staggered sites, and so changes S by ±1, so that 
H 1 } = 0. Thus H(-y) = (-l) s H(y)(-l) s . 

This observation means that m(y) is a polynomial in y 2 when appropriately normalized. Then 
in the expansion Q we have (—l) s \ty a ) = (— l) a |^ a ). States with even a are therefore orthogonal to 
those with odd a, and so m(y) can only depend on y 2 . For finite n, it is a finite polynomial; one needs 
to iterate only a finite number of times to determine the full wavefunction. The expectation value of 
any operator that commutes with (— l) s is therefore an even function of y, while any operator that 
anticommutcs is an odd function. 

Let us give a few explicit examples. For any state |a), define its coefficient in the ground states 
by ip a = ($ c \a) and Xa = (\£°|a). For N = 6 sites we find the components 

^ooo»o» — Xooo»o» — 

^•oo»oo 1 X»oo»oo 1 H~ 2y 

^•ooo»o ~~ y x»ooo»o y 

For N — 9 sites we have 

^ooo»o»o»o y Xoootototo ■ " y 

V'ooo^o^oo* — : ~y Xoooiotoo* — y (1 ~\~ 2y ) 

^•oo»oo«oo — - 1 ~h 4y X»oo»oo»oo — 1 ~h 2y -h 2y 

^.ooo.o.oo = -y(l + 3y 2 ) x»ooo.o.oo = -y(l + y 2 ) 

^•ooo»oo»o = y (l ~\~ 4y ) X* ooo * oo# ° = y 

^•oooo»o»o = ~2y^ X*oooo»o»o = 

V^o«oo«oo»o — : 2y (l -\- Ay ) Xo»oo»oo»o — ■ 
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The remaining ip a and \a are found by acting with T 3 and P and arranging the states into appropriate 
multiplets. Note that these coefficients are indeed even and odd in y depending on the value of 
(— 1) |a). With our sign convention P\0) = (— 1)™ +1 |0), P acting on all these states with fermion 
number n is found simply by inverting the picture; no extra minus signs occur. For example, 
P| • o o • o o) = ooiooi). Thus any state invariant under inversion does not exist for odd 
parity. 

3.3. Special properties of the ground state 

As discussed in the introduction, the ground states of these supersymmetric models possess a number 
of remarkable properties. We have above detailed some properties that can be derived using the 
supersymmetry. In this subsection we detail some special properties that do not automatically follow 
from the supersymmetry. They presumably are a consequence not only of the supersymmetry but the 
underlying integrability of the chain. However, as with the analogous results for the XYZ chain in 
[TU1 ITT1 [T2l IT5] , these properties are not derived but rather discovered by analyzing the explicit ground 
states. 

We define the polynomials characterizing the even- and odd-parity ground states with n fermions 
and N = 3n sites by 

m n (y 2 ) = (* c |^( ))/2, 
m_„(j/ 2 ) = (*°|*°(0)}/2 . 

These are normalized so that m„(0) = m_ n (0) = 1. Once these polynomials are known for a given n, 
the rest of the ground state follows simply by the interaction procedure. 

In fact, the following observation allows us to determine m± n iteratively on the computer using 
Maple, so that one does not need to simultaneously iterate while looking for m± n . Namely, with the 
normalization that m± n (0) — 1, we find (up to n = 8) the surprising fact that the norm of the ground 
states is also given simplify in terms of the m± n and rn±r n +iy. 

(* e (y)\*°(y)) = 2m„(y 2 )m_ (n+1) (y 2 ) , 

= 2m_„(y 2 ) m n+1 (y 2 ) . (7) 

This means that we can determine the polynomials for the ground state with n + 1 fermions simply 
by computing the norm of the ground states for n fermions. The first few polynomials are given by 

m 1 (y 2 ) = m 2 {y 2 ) = 1, m 3 (y 2 ) = 1 + 4y 2 , 

m A {y 2 ) = 1 + 5y 2 + 8y\ m 5 {y 2 ) = 1 + lly 2 + 42/ + 80y 6 + 64j/ 8 , 
m_i(y 2 )=l, m_ 2 (y 2 ) = 1 + 2y 2 , m_ 3 (y 2 ) = 1 + 2y 2 + 2y\ 
m_ 4 (y 2 ) = 1 + 7y 2 + 12y 4 + Uy 6 + 8y 8 , 
m_ 5 (y 2 ) = 1 + 9y 2 + 30y 4 + 46y e + 54y 8 + A2y 10 + 16y 12 . 

With the ground states in hand up to n = 8 (so that we have m± n (y ) up to n — 9), we have 
made the following observations: 

• All ip a and Xa are even or odd polynomials in y with integer coefficients. Note that the iteration 
procedure only requires that the coefficients be rational; the fact that they are integers is special. 

• All the integers in the polynomial for a given ip a or \a have the same sign. These signs are 
opposite for any two configurations related by moving a single fermion by one site. 

• The degree of m n {y 2 ) as a polynomial in y for n positive is \_(n — 1) 2 /2J, whereas the degree of 
m ~n{y 2 ) is [n 2 /2\, where [a\ is the largest integer less than or equal to a. 

• The polynomials m„ appear elsewhere in the ground state. Namely, the projection of the parity- 
even ground state onto the fully-polarized state (its ground state at y — > oo) is 

(V\V c (y)) = 2y n m n (y 2 ) (8) 

Moreover, the projection on states that are fully polarized except for one particle is —y n ^ 1 m n (y 2 ). 
This property can directly be confirmed by perturbation theory around the point y = oo. 
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• The projection of the parity-odd ground state onto the singlet state is 

(S\^°(y))=2m n+1 (y 2 ). 

• When a is a "fully squeezed" state like • o»o»oooo, \ a and ip a are monomials. 

• The polynomials m n (y 2 ) factorize into two polynomials in y 2 with integer coefficients when n is 
odd, while m_„(y 2 ) factorizes when n is even. For example, ms(y 2 ) = (1 + 2y 2 )(l + 9y 2 + 24y 4 + 
32y 6 ) and m-^y 2 ) = (1 + y 2 )(l + 6y 2 + 6y 4 + 8y 6 ). Below we will give a change of variables that 
will generalize the factorization to all n. 

• Evaluating the polynomials at the critical point y — 1 gives integers with interesting combinatorial 
properties [IB], very similar to what occurs in the XXZ chain at A = —1/2 (see e.g. [BJ). We 
detail some of these properties in the appendix. 



4. Duality and scaling of the order parameters 

In our earlier paper [5] we described how to find simple but exact expressions for several expectation 
values, and noted that there seemed to be a duality between large and small y, exact even on the 
lattice. In this section we review and expand upon these results. In particular, we find a useful change 
of variable that will allow us in the next section to find an explicit recursion relation for the m± n (y 2 ). 



4--1- Order parameters 

There are three distinct one-point expectation values in this supersymmetric chain. All three play 
the role of order parameters, distinguishing the weak-coupling (y < 1) phase from the strong-coupling 
{y > 1) phase. The staggered densities in the even- and odd-parity ground states are the expectation 
values of the operator S from ©: 

Using the derivation in section 13.21 these must be functions of y 2 . Our observation in section 13.31 
implies that these are ratios of polynomials in y 2 with integer coefficients. 

From the asymptotic ground states found in section 13.11 it is easy to see that 

^(0) = P°(0) = p°(y 2 = oo) = 0, £(oo) = l. 

Translation invariance at the critical point y = 1 constrains the expectation values, because here the 
combinations |\E' e ) ± l^ ) are eigenstates of the translation operator T. The eigenvalue of S in these 
eigenstates must be n/3, so this requires 

^(l) + p°(l) = 2/3. 

However, it does not follow that the difference of the two is zero at the critical point: S and T do not 
commute, so S has off-diagonal matrix elements between translation eigenstates. In fact, it not zero, 
as we will detail in section I4T31 below. 

In [5] we described how various expectation values were scale free. This means that the coefficients 
in the perturbative expansion around completely ordered points y = 0, oo are independent of N up to 
order y ±2n . Let us illustrate this here with the expansion of p e (y 2 ) + p°(y 2 ) around y 2 = oo using the 
exact ground states: 

28 V° + o(y- 12 ) 



P° 2 (y 2 ) 


+ p°2(y 2 ) 


= i 


pl(y 2 ) 


+ P3(y 2 ) 


= i 


Pl(y 2 ) 
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It is thus natural to conjecture that this pattern persists for all N, and so we can sum the series to 
give 

pUv 2 ) + pUv 2 )^y^ 2 , «>i. (9) 

The fact that this expression gives the correct value of 2/3 at the critical point y = 1 is a strong 
indication that this formula is both exact and valid for all y > 1. The expansion around y = has 
the same scale free property, and summing this series gives 

pUv 2 ) + pUv 2 ) = - - 2 = | , y < i (io) 

The two expressions are indeed continuous at y = 1, but the second derivatives are different, as one 
would expect at a critical point. 

Another order parameter involves the fermion densities on the sites 1, 3, 4, 6, . . . Because overall 
fermion number is conserved, it is easy to see that 

F = S + TST' 1 +T~ 1 ST 

is a constant. Therefore to get something new, we should look at the difference of the latter two 
operators. The difference is odd under parity, so it maps one ground state to the other. We define 

_ (^(TST-t-T-^m 

Rewriting this in terms of translation eigenstates at the critical point and exploiting the fact that the 
translation eigenvalues there are (— l)™ e ±I7I 7 3 gives 

rn(l) = ^ (/£(!)- /£(!)) • 

The order parameter t„ as well as the individual p° n are scale free as well, and summing their 
expansions as above gives conjectures analogous to ([9]) and (flO|) . We will give their expressions below, 
after finding a useful reparametrization of y. 

4-2. Duality 

Several considerations suggest that there is a duality between strong (y > 1) and weak (y < 1) 
coupling. Because the sine-Gordon field theory is the same for either sign in front of the off-critical 
perturbation, the scaling limit of this lattice model should be the same for y above and below the 
critical point y = 1. Moreover, as we will detail in section [5l the XYZ chain in the same universality 
class possesses an explicit duality symmetry. We will show in this subsection that there indeed is 
exact duality relation between certain combinations of the order parameter expectation values in our 
fermion chain. This duality relation is exact even for a finite number of sites. 

To exhibit this duality, it is very convenient to parametrize the staggering by a different variable. 
Examining the exact expression (0) suggests that we find a variable that removes the awkward square 
roots. An obvious choice is £ = yT+ 8y 2 , and indeed in the next section we will utilize this variable, as 
it makes the connection to the XYZ chain transparent. To understand the duality, £ is less convenient, 
because the critical self-dual point occurs at £ — 3. Thus in this section, we utilize the variable v 
defined by 

o 1 — V , . 

» " (TT^ ' <12) 

The critical point occurs at v = 0, with the y > 1 regime corresponding to — 1 < v < and the 
< y < 1 regime corresponding to 1 > v > 0. We will show that the duality transformation is simply 
v — > — v. 

To find which quantities are invariant under the duality transformation, let us first rewrite the 
asymptotic expressions for the order parameters in terms of v using (|12[) . By slight abuse of notation, 
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we will keep the same symbols, i.e. write p c n = p a n (v). The new parametrization allows us to combine 
([H]) and (fTU|) into a single expression valid in both regimes: 

2 8 v(v + 3) 

3 ~ 3(3+M) 2 

Obviously, this is not invariant under v — > —v. However, it does suggest that we study the quantity 

By analyzing the exact expectation values for n up to 8, we find that /3 n {v) is invariant under duality 
for finite values of n as well: 

Pn{v)=pn{-V) . 

Thus it is natural to conjecture that j3 n is a function of v 2 for all n. Of course, if one reverts to the 
original parameter y, it remains a function solely of y 2 as well. 

The other expectation values discussed above also have nice expressions in terms of v. By summing 
the scale-free series, we have 



Poo( u )-PooW = — 5\2 ' 0>U>-1 (14) 



y/v (v + 3){ V - 1) 

(«-3) 

while p%j(v) — p^o{v) = for 1 > v > 0. Combining this with (fT3)) gives for the product 

p c o>K>) = (^ 2 

for all v . This suggests we examine 

ln{v) = Pl(v)pl(v) . 

We have checked up to n — 8 that this quantity is indeed self-dual: j n (v) — 7„(— w). 
For the off-diagonal order parameter defined in (fTl~T) . we sum the series to get 

Too{v) = -^r, 0<V<1 (15) 

while Too = for > v > — 1. This has a square- root singularity at the critical point similar to (Ti"4"|) . 
and so suggests that we examine the product 

5 n (v) = (pl(v)-p° n (v))T n (v) . 

We have checked up to n = 8 that this quantity is indeed self-dual: S n (v) = S n (—v). 

4-3. Critical exponents and scaling functions 

We can extract various critical exponents from the exact expressions we have found. These match 
what is expected from conformal field theory, and so provide additional evidence that our conjectures 
for the N — > oo limit are correct. Moreover, they allow us to define appropriate off-critical scaling 
functions. A scaling function has a well-defined limit as N — > oo. To find one in an interacting 
system, the "bare" parameters (those in the Hamiltonian) must be rescaled. Here we will show that 
the appropriate rescaling is to plot quantities in terms of n 2 ^v. 

Let us start with the behavior precisely at the critical point. As we showed above, p^v = 
0) + Pn( v — 0) = 2/3 exactly for any n. Using the exact ground states, we find that for 
n = 2, 3, 4, 5, 6, 7, 8, p c n (v = 0) - p° (v = 0) is 

8 7 14 13 208 988 1976 



15' 15' 33' 33' 561' 2805' 5865 

A little trial and error shows that this sequence is given by 

e/ ^ o, n s 4 " 3fc - 2 4 r(2/3)r(l/3 + n) 
p n (v = 0)-p n (v = 0) = gll^ri = 3r(l/3)r(2/3 + n) 



k=l 




If we assume this applies for all n and then use Stirling's formula, we find that for large n 

-*(..B,. B -V»^( 1 -^ + ...). 

The fact that this vanishes as rT 1 ! 3 indicates that in conformal field theory this is the expectation 
value of an operator with dimension 1/3. There indeed is an operator of this dimension in the 
superconformal field theory identified in |3|4l[T4] as describing the scaling limit at the critical point. 

Going off the critical point, we see that both p%,{v) — p%o{v) and t^v) have square-root 
singularities in v as the critical point v = is approached. Since at the critical point, both quantities 
scale with system size as n^ 1 / 3 , the standard scaling argument indicates that the operator perturbing 
the theory away from the critical point has dimension 4/3. This argument is as follows. Perturbation 
theory around the critical point gives 

r n (v) rj A{n) (l + ax(n)v + a 2 (n)v 2 + . . .) . 

For n large we have A(n) oc 1/n 1 / 3 , and for v small we have limn^oo r n (y) oc \/v. The only 
way this is consistent for n large is for the coefficients ak(n) to diverge as n 2k / 3 . To ensure an 
appropriate scaling limit, we define a rescaled variable w = n 2 ^ 3 v that we keep finite as n — > oo. This 
variable w has scaling dimension 2/3, so if we consider an effective Hamiltonian near the critical point 
H = Hcrit + w J dx<f>(x), the operator $(a;) must have dimension 4/3. This indeed is the dimension 
of the only relevant supersymmetry-preserving deformation of this superconformal field theory |15j . 

Knowing these exponents, we can now find scaling functions that remain finite and converge 
quickly in the critical region as n — > oo. Instead of depending on n and v independently, they are 
functions of the single variable w — n 2 ^v. For example, we define 

g(w)= lim n 1/3 (T n (n- 2/3 w) - Tooin^^w)) . (16) 

To illustrate this, we plot (r n — t^) as a function of v on the left of figure[T]for sizes 9, 15, 21. Although 
these go to zero in the trivially solvable limits v = — 1 and v = 1, we see the effects of finite size strongly 
at the critical point v = 0. When we plot instead the rescaled version g(w) on the right side of figure 
[TJ we see an almost-instant collapse to a single curve. Thus we obtain a numerically-accurate scaling 
function from a very small system size. This strong collapse almost certainly is a consequence of the 
supersymmetry. 

5. Exact properties of the polynomials characterizing the ground state 

In this section we explain how the simple change of variables from y to £ = y/l + 8y 2 makes obvious 
a direct relation to polynomials that were studied by Bazhanov and Mangazeev [101 [HI H2], and 
Razumov and Stroganov |13j . in the framework of a particular XYZ chain. They turn out to be 
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related to classical integrable equations and the Painleve VI transcendent. This relation proves to be 
very fruitful, because several properties of the our chain can be deduced from already-known features 
of this XYZ chain. Here we discuss how the polynomials m± n (y 2 ) characterizing the ground state 
satisfy a recursion relation in n. Like the order parameters discussed above, they are scale free and 
obey a duality. The scale-free property in particular allows us to find an explicit expression for their 
n — > oo behavior. 



5.1. XYZ chain and the Painleve VI tau- function hierarchy 

Let us first recall some facts about the particular XYZ chain which possesses some striking properties 
[HI E3 ED H H3 GO • The XYZ Hamiltonian is 

1 L 

Hxyz = -- J^l+i + J v°>t+i + J*°k°i+i > ( 17 ) 
fc=i 

where the Pauli matrices act on a the two dimensional Hilbert space at each site. As above, 
we take periodic boundary conditions. The special properties occur along the line of couplings 
JxJy + JxJz + JyJz = 0. A convenient parametrization of the coupling constants along this line 
is given by 

Jx = 1 + C, Jy = 1 - Jz = (C 2 - l)/2 • (18) 

This model is an off-critical deformation of the XXZ model at A = J z /J x = — 1/2 (the C = case). 
In the conventions of our earlier paper [5], £ = 3/s. 

An obvious symmetry of the XYZ spectrum comes by permuting the Pauli matrices. Along 
the line (|T5|). this results in two useful transformations. A rotation of all spins by tt/2 around 
the z-axis exchangs the a x a x - and <7 y cr y -terms in Hxyz- This is equivalent to the transformation 
( — > — £. Conversely, rotating all spins by tt/2 around the y-axis exchanges the a x a x - and er z <7 2 -terms. 
Equivalently this can be done by the homographic transformation 

C^§±f <») 

and rescaling energy by (C — l) 2 /4. Both these transformations are symmetries of the spectrum of 
Hxyz, up to an overall rescaling. 

When the number of sites L is odd, all the eigenvalues of Hxyz are doubly degenerate under 
spin reversal. As Hxyz commutes with both the operators I = J\k=i °fc anc ^ ^ ne spin-reversal 
operator R = Ylk=i ^fe one ma y decompose the corresponding subspaces into eigenfunctions with 
/|$±(0) = ±|$±(£)). These are related by spin-reversal -R| ( I > ±(C)) = I^t(C)): an d thus one may 
concentrate on one of them, say |$_(£)). 

Let us now turn to the ground states of the model. For an odd number of sites L = 2n + 1 the 
ground states |$±(C)) have energy E = —L(C, 2 + 3)/4 exactly [TTJ [7]. Like for our fermion chain, 
consider a decomposition of the wave function into spin-configuration states | <&_(£)) = J2 a <f) a (C) | a) . 
Given that the Hamiltonian itself and the ground state energy are quadratic polynomials in £, 
one can choose the normalization in such a way that all components are polynomials in £. An 
overall constant normalization can be fixed by demanding that the zero-order term for configurations 
a =tt ' ' ' til ' ' ' I containing a block of n (resp. n + 1) spins up for n even (resp. odd) equals 
one. Then it turns out that all polynomials 4> a (C) have integer coefficients and definite parity. Some 
of them enjoy special properties: the component in front of the fully-polarized state \4- ' * ' X w ith all 
spins down is given by 

hi-i(0 = C" ( " +1)/2 sn(C 2 ) (20) 

where the s n (z) are polynomials in z of degree [?i 2 /4j that solve the non-linear recursion relation 
2z(z - l)(9z - l) 2 (lns„(z))" + 2(3z - l) 2 (9z - l)(lns„(z))' 

+8 ( 2 „ + 1)2 £n+i(£K-i(£) _ (4(3n + 1)(3n + 2] + (Qz _ 1)n(5n + 3) ) = o (2i) 

s„{zy 
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with initial conditions sq(z) = Sx(z) = 1. In fact, these are a special case for the Hirota equations 
for the r-function hierarchy corresponding to Picard solutions of the Painleve VI differential equation 
[TBHH]. Moreover, these polynomials also appear in the computation of the Q-function for the ground 
state of the corresponding eight- vertex model. 

There is a number of remarkable properties of the XYZ chain involving the polynomials s n (z); for 
an extensive list see jT2l [13] . As an example let us consider the norm of the wavefunction. Finite-size 
computations show that it is given in terms of the polynomials s n {z) by 

||$_(C)|| 2 = (4/3)"C n( " +1) Sn(C 2 )s-n-i(C 2 )- (22) 

Using the invariance of the ground state sector under (fTTJf it is not difficult to show that (|22p is 
invariant under the homographic transformation modulo rescaling by a ^-dependent function. After 
a little algebra, one finds the invariant 

S n (C 2 )S-n-l(C 2 )=(^) nin+1) Sn\ j^S |.s_„_, I 1^-1) I . I 23 I 





2C 

5.2. The polynomials m± n (y 2 ) and Painleve VI 

We noted in section 13.31 that m n (y 2 ) factorizes into a product of polynomials in y 2 with integer 
coefficients when n is odd, and m_ n (y 2 ) factorizes when n is even. However, as we detail in the 
appendix, they factorize for all n at the critical point y = 1. Thus one might hope that they factorize 
for all y with an appropriate change of variable, and indeed changing to £ or v allows this factorization. 

This change of variables allows us to find a much more dramatic result: the m± n (y 2 ) are simply 
related to the analogous polynomials s n (z) in the XYZ chain! By examining their explicit form up 
to n = 8 in terms the new variable £, it is readily apparent that they are related to the s n (z) by the 
following equations: 

/-L(«-1) 2 /2J s ,(Y-2\ 

/d« 2 /2j s 

^ 2 -(C 2 -D/8)= 2( i lKB+2) /-y . (25) 

with s„(0) = 1 for n > 0, and s„(0) = (3/4)-™- 1 for n < 0. 

The most important consequence of this identification is that the polynomials m± n (y 2 ) can be 
obtained from a recursion relation. Moreover, the ground states in the fermion and XYZ chains have 
analogous special properties: for example, (|2"U)l is analogous to ([5]). 

We can use the knowledge about properties of the s n (z) in order to obtain ground state properties 
of our fermion chain. When the transformation (jTHJ) is rewritten in terms of v — (3 — £)/(£+ 1), this 
is precisely the v — > — v duality discussed in section (I4.2[) . The square norm of the ground state 
wavefunctions of the fermion chain, given in ([?])■ are quite similar, although not identical to (|22p . 
However considering the combination of (0, (|24|) and (|25|) . it is a straightforward consequence of 
(|23ll that the product of the ground state square norms N n (v) = \\^ c (y)\\ 2 \\^ (y)\\ 2 has the following 
invariance property 

/ I 1 \ n(n+l) 

Nn(~v)= *n(v)- (26) 



Given this invariance property it is not surprising that all quantities studied in section 14.21 with nice 
duality properties involve both parity sectors. 



5.3. Scale-free properties 

Knowing that the m± m (y 2 ) are related to the Painleve VI equation teaches us a great deal about the 
ground states of our fermion model. Here we show how the lessons we have learned in the analysis 
of the fermion model give new insight into the polynomials s n (x). In particular, we show how they 
obey the scale- free property discussed in [5] and above in section 14.21 Combining this with the Hirota 
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equations ([2"Tj) allows us to derive the scale-free coefficients explicitly. We can then sum this series by 
approximating the Hirota equation, and so find an asymptotic expression for the s n themselves. 

To this end, let us introduce another reparametrization x = (1 — £)/(3+C)- Then define a rescaled 
version of s n via 



S4x) = s ,4l/C 2 )(l+x)^(l-:r)^(l + 3x)^/ S „(l) . (27) 

where s„(l) = 2™("- 1 )/ 2 for n > 0, and s„(l) = 3-«-i2(™+ 2 )(™+i)/2 f or n < 0. The Taylor series 
expansion of S n (x) around the trivially solvable point x — is scale free: for n > the first n terms 
in this expansion are independent of n, while for n < the first — n — 1 terms in this expansion are 
independent of n. The scale-free part of the expansion is 

S{x) = 1 + x + 2a; 3 - x 4 + 9x 5 - 6x 6 + 54x 7 - 36x 8 + ... 

We can determine all the coefficients and a summed form directly from (l2~Tj) : the S n (x)'s are 
solution to the recurrence equation 

p(x)0nS n (x))" + q(x)(ln S n (x))> + 2(2n + I) 2 (^^±lM^iM _ ^ = r{x)> 

with the coefficients 

2x(x-l){x + l) 2 (3x-l) 2 



p(x) = 
q{x) 



l + Zx 

(x + l)(3x - 1)(1 + 8x - 22x 2 + 45x 4 ) 



(l + 3x) 2 
r{x) = -3x 3 - x 2 + 3x - 1. 

Let us now make the assumption that the limit Sodx) = linin^oo S n (x) exists and that we can neglect 
the term oc (2n + l) 2 (i.e. the ratio S n+ i(x)S n -i(x) / S n (x) 2 — 1 decays faster than 1/n 2 , at least for 
small x). Thus we must solve the differential equation 

pix^lnSooix))" +q(x)(\nS QC (x))' = r{x) (28) 

with initial conditions Soo(x = 0) = 1, S' oc (x = 0) = 1. The unique solution is given by the simple 
function 

(1 + 3x) 5 / 24 

SOO{X) = (l-. T )3/8(l + . T )l/8 (1 „ 3a , ) l/24 ( 29 ) 

which indeed generates all the scale-free coefficients. 

An illustration of this function compared to S n (x) with n = ±50 is given in figure Ufa). The 
recurrence relation allows us to compute the first corrections to the scale-free form as well. We have 
S n (x)/Soo(x) = 1 + b n x n+1 + c n x n+2 + ... Upon inserting this into the equation for S n (x), and 
expanding around x = we find the recursion relations (n + l)b n = 2(2n + l)6 n _i, bo — —1 and 
(n + 2)(2n + 3)c„ — 2(2rt + l) 2 c„_i — 4b n = 0, do = 1/2 which are readily solved: 

, 1 / 2n + 2\ 1 / 2n + 2 \ . . 

K = —2\ n + 1 J' Cn= ^T2{ n + 2 J' ( 30 ) 

Hence the first correction is a binomial - a quite common feature that we observed for other quantities 
as well. 

The preceding procedure relied on an expansion around the trivially solvable point x = 0. As we 
did with the order parameters, we can also understand the scaling around the critical point x = 1/3. 
From the asymptotic formula (|29[) . we see that Soo{x) displays a power-law singularity with exponent 
— 1/24. Conversely, for finite n we find at this point that 

gn+i(l/3) _ 2 4 ("+ 1 ) S » +1 (l/9) S »(l) _ 2 4 »+ 3 n!(3n + 2)! 
S„(l/3) 3 2 »+5/2 s „(i/9) s „ +1 (l) 3 3 ™+5/2((2n+l)!)2 [ > 

The value s n (l/9) necessary here can directly be inferred from the Hirota equations by setting 
z = 1/9. Upon application of Stirling's formula for large n, it is not difficult to derive from this 
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expression the asymptotic behavior S n (x = 1/3) ~ n 1 / 36 for large positive n. For large negative 
n the scaling is identical because of 5_„_i(l/3) = s/3S n (l/3), Thus to obtain something finite in 
the n-)oowe obtain again the result from section |4~31 that near the critical point, the deformation 
parameter (here x — 1/3) must be rescaled by a power of n 2 / 3 . In the context of differential equations, 
such a relation between properties in the two limits is known as a "connection formula" . The 
rescaled functions a n {t) — 7Z. _:L ' 36 5 n (l/3 — n~ 2 / 3 t) are illustrated in figure HJb). As before, we 
observe a remarkable collapse for large n, this time however without subtraction of any infinite- 
volume term like n~ 1 / 36 5 00 (l/3 — ri~ 2 / 3 t). It would be very interesting to find the scaling functions 
cr±(t) = lining-toe a n (t). A detailed analysis of the Hirota equations [23] yields at least the first two 
Taylor coefficients 

CT±(<) = C± ( 1± 4TW|^ ^)' (32) 

where G+ = G(3/2) 2 /(G(4/3)G(5/3)) = 1.013993 . . . and G_ = V3C+. Here G{z) denotes the Barnes 
G-function, defined through the functional equation G(z + 1) = T(z)G(z). 

6. Conclusion 

We have studied an off-critical staggered model for supersymmetric lattice fermions with hard-core 
exclusion, and found some striking properties of its ground state. We showed how the supersymmctry 
allows the entire ground state to be characterized by a single polynomial, and related this to solutions 
of the Hirota equations, which earlier arose in an XYZ chain in the same universality class. We 
saw how a variety of special properties allowed the determination of exact expressions for the order 
parameters and critical exponents. 

The precise connection between the supersymmetric fermion model and the XYZ chain is 
surprising for several reasons. First, it is known that the critical point of the fermion chain can 
be transformed to the XXZ chain at A = —1/2 with twisted boundary conditions and an even number 
of sites [U H]. Even though there is a relation between the cases L = 2n with twisted boundary 
conditions and L — 2n + 1 with periodic boundary conditions for this XXZ chain |20j , what happens 
off criticality remains mysterious. In particular, no XYZ chain with twisted boundary conditions that 
remains integrable is known, since there it has no obvious U(l) conservation law. Note, however, the 
fermion chain does possess a t/(l) fermion number conservation law. Another distinction between 
the two is in the dualities. For the XYZ chain, the duality transformation (fT9|) is related to simple 
unitary transformations on the Hilbert space (spin rotations) and obvious invariance properties of 
the corresponding Hamiltonian. In the fermion case however, the transformation on the coupling 
constant y is highly non-linear, and a realization of the duality on the Hilbert space through a unitary 
transformation seems to be less evident. 
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These observations give a strong hint that there remains much structure to be uncovered in both 
the fermion and XYZ chains. As the latter are related to the eight-vertex model at the special crossing 
parameter r) = tt/3, powerful theoretical tools such as Baxter's famous TQ-equation can be used for 
their analysis. For example, it allows the exact evaluation of the nearest-neighbor spin correlators in 
finite size |23] . In the case of the fermion chain, similar relations remain to be discovered. 
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Appendix A. Combinatorial properties at the critical point 

Given the connection between the critical model and the XXZ chain at A = —1/2, it is natural 
to expect that the values rrv^(y = 1) are related to combinatorial numbers enumerating alternating 
sign matrices and related objects. By explicitly finding the ground states at y = 1, some of these 
combinatorial properties were explored in [16] . Here we analyze the particular case of N — 3n and 
periodic boundary conditions in more depth. 
First, we notice that at y = 1 we find 

m_ 2 fc(l) = 3m 2 fe(l) and mi_ 2fe (l) = m 2fe _i(l), 

for k = 1,2,.... Therefore it will be sufficient to concentrate on the parity-even sequence m n (l). For 
odd n = 2k — 1 we thus find m 2 fc_i(l) = 1, 5, 198, 63206, 163170556, . . . This sequence corresponds to 
the numbers Ajjjj (Ak) of 4fc x 4fc alternating sign matrices of type UU [3T] . It is known to factorize 
into 



where 



m 2fc _i(l) - Auu(4k) = A v (2k + 1) x A uv (4k; 1, 1, 1) (A.l) 



A foirMU 1 TT (6.7 -2)!(2. 7 -l)! 

Ay(2fc + 1) =^11 



«l = l(4i-l)!(4j-2)!' 



2/c 2A;-|"-1 q ■ Q * i o 

A TO (4fc;l,l,l) = (-3) fe2 2 2fe n II 



j - i + 2n + 1 

Here Ay{2k + 1) is the number of alternating sign (2k + 1) x (2k + 1) matrices symmetric about the 
vertical axis, also the number of (2k) x (2k) off-diagonally symmetric alternating sign matrices, and 
Auu(Ak; 1, 1, 1) expansion coefficients of a generating function appearing in [2T | 122 ] . 

For even n = 2k we find the values m 2 fe(l) = 1, 14, 1573, 1427660, ... As is the case of an odd 
number of particles, this sequence can be factorized according to 

m 2k ( y = 1) = N 8 (k) x ^ ( ( 2 fc fc ~- 1 l) ' fc = 1 < 2 '--- (A- 2 ) 
with the numbers of cyclically symmetric transpose complementary plane partitions N$(k) given by 

m 

Mm tfftai^-ilL 
W = n (3j + 1) (4i)!(4i + l)!' 

and the well-known numbers ofnxn alternating sign matrices 

We emphasized the factorizations because it turns out that the subfactors of m 2 fc_i(y 2 ) and 
m -2k(y 2 ) are in one-to-one correspondence with the factors in (|A. 1|) and (|A.2[) when setting y = 1. 
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